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Abstract 

In this paper we introduce digital nets over finite abelian groups which con- 
tain digital nets over finite fields and certain rings as a special case. We prove 
a MacWilliams type identity for such digital nets. This identity can be used to 
compute the strict i-value of a digital net over finite abelian groups. If the digital 
net has iV points in the s dimensional unit cube [0, l] s , then the t- value can be 
computed in O(iVslogiV) operations and the weight enumerator polynomial can 
be computed in 0(Ns(log N) 2 ) operations, where operations mean arithmetic of 
integers. By precomputing some values the number of operations of computing the 
weight enumerator polynomial can be reduced further. 



1 Introduction 

Digital nets are point sets {xq, . . . , cc^-i} in the s dimensional unit cube [0, l) s with 
N = b m points (where N,b > 2 and m > 1 are integers), whose construction is based on 
linear algebra over a finite field (or more generally finite rings) [H QUI US] . The aim of 
these constructions is to obtain highly uniformly distributed point sets. Such point sets 
are useful as quadrature points for quasi-Monte Carlo rules 

1 N ^ 
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which are used to approximate integrals J, Q 1 - )S f(x) dx. The Koksma-Hlawka inequality 
states that the error is bounded by the variation of the integrand / times the discrepancy 
of the quadrature points [U Chapter 2]. One measure of the distribution properties of 
digital nets is the concept of (t, m, s)-nets in base b [U [TS]. A point set P in [0, l) s 
consisting of b m points is a (t, m, s)-net in base b if every interval of the form 

s 

n 

i=l 

for every choice of integers < a, t < b di and integers d±, . . . , d s > with d\+- ■ -+d s = m—t, 
contains b l points. Thus the t- value measures the quality of (£, m, s)-nets (where smaller 
is better) [2J [15]. If t is the smallest integer such that the above property holds, then P 
is a strict digital (t, m, s)-net and we call t the exact t-value of P. 

Many explicit constructions of (t, m, s)-nets with good quality parameter t are known, 
see [HI EH El [25] or [3J Chapter 8]. These are based on the concept of digital nets which 
we recall in the following. For a positive integer b, Z 6 = {0, 1, . . . , b — 1} denotes the 
residue ring modulo b. For a prime power q, ¥ q denotes the g-element finite field. Let 
Ci, . . . , C s G Fg Xm benxm matrices over a finite field ¥ q (with n > m). Let <p : Z g — > ¥ q 
be a bijection. For < / < q m let 

m— 1 

/ = l r q r , where l r G Z g , 

r=0 

be the base q expansion of I. Let / = (f(lo), ■ ■ ■ , f(l m ~i)) T G and 
where y jtl = (yj,i,i, Uj,i,n) T G F™. Then we define 

, , l P~ 1 (yj,l,n) f i / ■ / 

£z j = — h 1 — — , for 1 < i < s. 

q q n 

The point set X\ = (2^1, . . . , for < I < q m is called a digital net (over F ? ). A digital 
net which is a (strict) (t, m, s)-net is called a (strict) digital (t, m, s)-net. 

Let Cj = (cT l5 . . . , cJ m ) T , i.e., Cj^ denotes the kth row of Cj. Then the condition that 
Ci, . . . ,C S generate a digital (t, m, s)-net over F g is equivalent to the condition that for 
all nonnegative integers d±, . . . , d s with d\ + • • • + d s = m — t, the vectors 

c l,l) • • • ) C l,dl> • • • ) C s,li ■ ■ ■ i c s,ds 

are linearly independent [H Theorem 4.52]. Using this definition it is expensive to 
compute the exact t- value (i.e., the smallest value of t for which the linear independence 
condition holds) since in general many linear independence conditions need to be verified. 
A direct computation of the t-value based on the linear independence properties of the 
generating matrices over a finite field is presented in [18]. In this paper we show how the 
exact t-value can be computed without checking any linear independence condition. This 
is done by using a Fourier inversion method which yields a fast algorithm for computing 
the exact i-value. Apart from the explicit constructions of digital nets, computer search 
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algorithms of (t, m, s)-nets are also useful in that they often yield digital nets with very 
small t- value PJJ [HI [21] . The usefulness of such search algorithms is limited by the size 
of the search space of digital nets and the computational cost of computing the t- value of 
some given point set [T8] . 

In this paper we generalize the concept of digital nets over finite fields jU [TH] or finite 
rings [10] to digital nets over a finite abelian group G. We also generalize the concept of 
duality theory of nets to finite abelian groups as studied in [U] for finite fields. In our 
context, the dual net is now defined via the dual group of characters of the finite abelian 
group G. We show how the quality parameter t of a digital net is related to the NRT 
(Niederreiter [2] and Rosenbloom-Tsfasman [2U]) weight in the dual group of characters. 

Further we introduce an algorithm which allows one to compute the weight enumerator 
polynomial H21 [271 EE] of the dual of a given digital net in 0(Ns(logN) 2 ) operations 
(see Algorithm [lj and the quality parameter t of a given digital net in 0(Ns log N) op- 
erations (see Algorithm [2]), where operations always mean arithmetic of integers. By 
precomputing some values or using a faster polynomial multiplication algorithm this con- 
struction cost can (theoretically) be reduced further. For instance, if one wants to compute 
the t-value of many digital nets (as in a computer search algorithm which generates many 
digital nets and then chooses the best one, as in jT8j 121] . or one wishes for instance to 
optimize the direction numbers of a Sobol' sequence as in [9]), it can be beneficial to store 
some values which might have to be computed repeatedly otherwise. It is also possible to 
balance the computational cost and storage cost in different ways, as will be clear from 
the result below. 

The main idea for the fast computation of the t-value is to interpret the worst-case 
error of integration in the Walsh space [3] in a different form. This yields another way of 
measuring the quality of a digital net (which is not based on the t- value). In this case a 
possible criterion is of the form 

1=0 fce{0,...,6 m -l} s 

where fi(k) is the NRT weight which measures the magnitude of k in some sense and 
walfc is the fcth Walsh function (for details see Section [2]). In the following we replace b~ l 
by a variable z, thus we consider 

- 1 + ^E E ^wal^O- 

(=0 fcG{0,...,6 m -l} s 

We view this expression as a polynomial in the variable z, i.e., we do not substitute any 
value for z, but rather aim at expressing this polynomial in the form 

sm 
a=l 

This polynomial is shown to be the weight enumerator polynomial of the dual net £>, see 
[H Definition 4.76] or Section [2] for details, namely we have a MacWilliams type identity 

fe m — 1 ms 

-Y, E ^ (fc) wai fc (^) = E^- 

1=0 fce{0,...,f> m -l} s a=0 
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Different but similar identities have previously been studied in coding theory, association 
schemes, orthogonal arrays, digital nets and so on [3 Q21 E7J [28]. In particular, a similar 
approach was considered in [27], but with a different aim. 

Here the aim is to compute the coefficients iV a for 1 < a < m. We show that these 
coefficients can be computed in 0(Ns(\ogN)) operations. On the other hand, since these 
coefficients are related to the t- value of the digital net, we can also compute the t- value 
of a digital net in 0(Ns(log N)) operations (see Algorithm [1]). 

A second result uses a Mac Williams identity in the reverse direction (see Algorithm |2]). 
This simplifies the computation of the t- value, but it does not yield the weight enumerator 
polynomial. The construction cost of this algorithm is 0(Ns log N) operations. 

We present some background in the following section. Section [2] introduces the notion 
of (t, m, s)-nets over finite abelian groups, and deals with the proof of the main results 
and gives some discussion of implementing the algorithm. We also include some results if 
the point set is not a digital net. In Section [3j as a proof of concept, we use our methods to 
compute the exact t-values of Sobol's sequence [2B] (as implemented in Matlab 2011a) for 
dimensions 3 < s < 22 and number of points 2 m with 2 < m < 25. In our experiments, 
Algorithm [2] is slightly faster than Algorithm [1] for computing the t- value (we did not use 
precomputation for both algorithms). In Section H] we introduce a further generalization 
of (t, m, s)-nets and digital nets, over finite abelian groups. We also study for which finite 
rings the ring-theoretic dual net coincides with the character-theoretic dual net used for 
digital nets over finite abelian groups. 

2 Digital nets and Walsh functions over finite abelian 
groups 

Let (G,+) be a finite abelian group of b elements. Let T = {z G C : \z\ = 1} be the 
multiplicative group of complex numbers of absolute value one. 

Definition 1. Let G be a finite abelian group. The dual group of G is defined by G* := 
Hom(G, T), namely the set of group homomorphisms from G to T, often called characters 
ofG. 

The following results are well-known [2U Part I]. Let G, G\, G2 be finite abelian groups. 

Lemma 1. 1. Let e be the exponent of G, namely, the maximum order of elements in 
G. Then, the image 0/ Hom(G, T) is exactly the cyclic group < ( e > of order e, 
where ( e := exp(27ri/e). 

2. G* is isomorphic to G as a finite abelian group, but there is no canonical choice of 
the isomorphism. 

3. A group homomorphism f : G\ — > G2 induces f* : G2 — > G\ by composition k G 
G* 2 ^ kof G G\. 

4- G — > (G*)* , x 1 — y (x** : k i—)- k{x)) gives an isomorphism, through which we identify 

g = (G*y. 

5. A morphism f is surjective (respectively injective) if and only if f* is injective 
(respectively surjective). 
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6. The dual of G\ x G 2 is canonically isomorphic to G\x G\. 

7. For any k G G* , the sum ^2 x&G k(x) is if k ^ 0, and if k = 0. Dually, for 
any x E G, the sum J2keG* M x ) ^ s if x ^ 0, an d if x = 0. 

Remark 1. T/ie notion of dual G* generalizes to locally compact Hausdorff abelian groups 
G by defining G* := Hom cont (G,T) (namely, the continuity is required for the characters), 
and the above lemma holds for this wider class (with some small adjustments), where finite 
abelian groups have discrete topology. See IW$ or [TI\ Chapter VII] for the basic of such 
groups and Fourier transformation on them. 

We return to our finite abelian group G. Define the direct product of denumerably 
many copies of (G, +) and the direct sum of denumerably many copies of (G*, +) by: 

X := G N = {x=(£i,6,---):&eG}and 

K := {k= («!,«2,...) G (G*) N : = for almost all i}. 

We remark that K and X are dual to each other, through the pairing •: 

00 

KxX^T, (ki,K2,...)«(£i,6,...) :=n«*(6) G T " 

i=i 

For a positive integer n, we have a truncated version 

n 

K n : = (G*) n , X n := G n , ( Kl , ...,K n )* (Ci, . . . e„) := J] G T. 

i=l 

These finite groups are again dual to each other. 

We fix bijections tp : Z b -> G and <j> : Z 6 -»■ G* with y?(0) = and 0(0) = 0. The 
bijection gives a bijection denoted by the same symbol by an abuse of notation: 

<p : N ->• K 

in the following way: We identify elements in IK with a nonnegative integer k G No 
with fe-adic expansion k — Ki + K 2 b + ■ ■ ■ + n m b m ~ l with G Then we define 
(j)(k) := (0(ki), . . . , 0, . . .) G K, which is a bijection. By restriction, we have a 

truncated version 

0:{O,l,...,6 m -l}^K m . 
Analogously we use the same symbol ip to denote 

cp : [0,1) -»X 

defined as follows. For a number x G [0, 1) with base expansion x = + ^fr -2 + ■ ■ • j 
we define <p(x) := (y?(£i), ^(£2)7 • • •) £ X. Here we use the finite expansion of x if x is a 
6-adic rational. Because of this reason </? : [0, 1) — >■ X is not surjective. However, both 
have natural measures (namely, Lebesgue measure and Haar measure) and tp preserves the 
measure, through which [0, 1) and X are isomorphic as a measured space when a measure- 
zero subset is removed from X. There is also a continuous mapping p~ l : X — > [0, 1] given 
by (£1,^2, • • •) ^ <£> _1 (£i)fr 1 + t f~ 1 {i2)^~ 2 + ' ' ' > which is surjective (and inverse to <p if 
measure zero sets are neglected). 
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The Walsh functions are then defined the following way: for k = ki + K2&+- • ■ + K m b m 1 
and x = + ^2^~ 2 + • ■ ■ we have 



G-walfcfx) = <p(k) • <f(x) = Jj0(«i)(<^(^))- 



i=i 



For the s- dimensional case, we again use the same letters 

0:f%-+K s , : [0, 1) S ^X S . 
For k = (ki, . . . k s ) and x = (x\, . . . , x s ) we define 



G wal fc (a;) = JJcwalfc^Xj) = (p(k) • tp{x) 



i=l 



In the following we recall the definition of (t, m, s)-nets in base b [U 

Definition 2. Let b > 2 and s > 1 be integers. A point set P = {x , . . . , a?( ) m_ 1 } c [0, 1^ 
is called a (t, m, s)-net in base b if for all nonnegative integers d\, . . . ,d s with di + - ■ -+d s - 
m — t the elementary intervals 

a,- a,- + 1 



n 



2=1 

Ut 



contain exactly b points for all choices of < a« < b 1 for 1 < % < s. 

Ift is the smallest value such that P is a (t, m, s)-net, then we call P a strict (t, m, s)- 
net and t the exact quality parameter (or t-value). 

Because of this definition, we may replace the coordinates of the points to be m- 
digit 6-adic rationals. We fix an n > m, and we consider point sets with coordinates 
being all n-digit 6-adic rationals. Namely, we consider only the image of tp^ 1 : X* — > 
[0, l) s . (For estimating the computational complexity we use n = m to avoid unnecessary 
complications.) 

We identify X* = (G n ) s with the set of s x n matrices (M s n (G), +), and similarly 
identify K s n = ((G*) n ) s = (M SjU (G*),+). For X E X* we write X = (&,,•) i<i<< and 

similarly K = (K itj ) i<j< s . 

l<j<n 

Definition 3. A subset P C X^ is a combinatorial (t,m, s)-net in base b if <£> -1 (P) is a 
(t,m, s)-net in base b. 

Note that this definition does not use the group structure of G. Note also that to 
define (t, m, s)-nets it is enough to consider the case n = m, but for the convenience in 
future use we choose to have general n. 

Definition 4. Let P C X* be a sub abelian group. We define P 1 - C as the kernel of 
(K n ) s = (X s n )* -> P*. In other words, 

P ± ■= {K E K s n I K • X = 1 for any X G P}. 

We call P 1 - the character theoretic dual net of P. 
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Let Y be a subset of {1, 2, . . . , s} x {1,2,..., n}. Let X* (Y) be the direct product of 
#(Y) copies of G with index in Y, namely: 

K(Y) :=m,Mh3)eY,^eG}. 

We have an obvious projection 

Pry = K ~> K( Y ), (kikj ^ (&i)(i,j)eY- 

Let P be a subset of X* , with inclusion map i : P — > X* . A function / : jSi — > S2 
is said to be uniform if #(/~ 1 (s2)) is constant for any s 2 G 5*2. In particular, it implies 
surjectivity if Si is non-empty. Conversely, if Si and 5*2 are finite groups and / is a 
group homomorphism, then the surjectivity of / implies uniformity, because any element 
si in the set f' 1 (s 2 ) gives a bijection / _1 (0) — > f~ 1 (s2) obtained by adding si. 

The point set P is Y -uniform if pr y o 1 : P — y X* (Y) is uniform. 

Lemma 2. Fix a positive integer m. A finite subset P C X^ 0/ cardinality b m is a 
(combinatorial) (t, m, s)-net if and only if the following conditions are satisfied: 
For any integers di, . . . , d s > 0, letY := Y(di, . . . , d s ) denote the union of 

{(l,l),(l,2),...,(l,rf 1 )},{(2,l),(2,2),...,(2, C / 2 )},...,{(s,l),(s,2),...,( S ,4)}. 

Then, for any di + d 2 + ■ • • + d s < m — t, the composition ?4X*^ "K^(Y) ^ s uniform, 
namely, P is Y -uniform. 

Proof. This follows from the definition of (t, m, s)-nets. □ 

The surjective homomorphism X^ — > X*(Y) induces an injective homomorphism 
(X£(Y))* -> (X*)*. Through the identification (X*)* = K s n = M SjTl (G*), (X£(Y))* C K s n 
is identified with 

(K(Y)T = {K= (Kv^Kij = for all £Y}c M s , n (G*). 

The following easy lemma links the (t,m, s)-net property of P to the minimum 
Niederreiter-Rosenbloom-Tsfasman (NRT) weight (see, for example, [U Section 7.1]) of 
P ± . 

Lemma 3. Let P C X^ be a sub group. For any subset Yd {1, 2, . . . , s} x {1, 2, ... , n}, 
the composition P C X* — > X* (Y) is surjective (namely P is Y -uniform) if and only if 
P 1 - Pi (X*(Y))* = {0}. In other words, if and only if any K G P ± with Ki j = for all 
(i,j)£Y isO. 

Proof. By m Lemma [T], the surjectivity is equivalent to the injectivity of X* (Y)* — > P*, 
which is equivalent to X^(Y)* D Ker(K* — > P*) = {0}, and the result follows from the 
definition of P L = Ker(K s n -> P*). □ 

For k = (aci, . . . , K n ) G (G*) n = K n which is not 0, we define 

n(k) = max{j|l < j < n, Kj ^ 0} 
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and for k = we define /x(0) := 0. For any K e (X^)* = M s , n (G*), let ki be the i-th row 
of K, and define 

s 

Ai(iO=X)A*(*i), (!) 

i=l 

which is nothing but the definition of the NRT- weight of K. For a subgroup Q C , we 
define its minimum NRT-weight 

minNRT(Q) := mm{fi(K)\K eQ,K ^ 0}. (2) 

If Q — {0}, we define minNRT(Q) = ns + 1. A direct generalization of the duality 
theorem for digital nets in [16] is the following: 

Theorem 1. Let P C X^ be a subgroup of cardinality b m . Then, P is a (t,m, s)-net if 
and only if minNRTlP- 1 ) > m — t + 1. In other words, m + 1 — minNRTlP- 1 ) gives the 
exact value oft for P, i.e., P is a strict (t,m, s)-net. 

Proof. By Lemma [2, the (t, m, s) condition is equivalent to the surjectivity of pr y o i 
(since uniformness is equivalent to the surjectivity for group homomorphisms) , and by 
Lemma [31 the condition is equivalent to that P 1 - n X* (Y(di, . . . , d s ))* = {0} holds for 
every d\ + ■ • • + d s < m — t, which is equivalent to that any K ^ with [J>(K) < m — t 
does not belong to P x , hence the condition of minimum weight of P^. □ 

2.1 A MacWilliams identity for digital nets over abelian groups 

In the following we define the weight enumerator polynomial. 

Definition 5. We define the weight enumerator polynomial WP P ±(z) of P with respect 
to the NRT-weight, which is a non-negative integer coefficient polynomial in variable z, 
by 

ns 

WP P i( 2 ) = zKK) = ^N a z a . 
This gives definitions of N a , i.e., N a is the number of matrices in P with NRT-weight 

a: 

N a = #{K EP ± HK) = a}. 

Thus we have 

minNRT(P ± ) = min{a|l < a < ns, N a ^ 0} 

and if N\ — ■ ■ ■ — N ns = we have minNRT(P J -) = ns + 1 (this means P 1 - = {0}). 
We introduce the notion of dual NRT-weight. 

Definition 6. We shall define the dual NRT-weight /i* on X n : for x = (£i, . . . , £ n ) £ 
X n = G n , we define jj*(x) as the index of the minimum nonzero coordinate, i.e., if x ^ 
we set 

H*{x) = min{i|^ ^ 0}. 

If x = 0, then (A*(x) := 0. 
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It is natural to define fi*(X) for X — (x\, . . . , x s ) G X^ by 

s 

i=i 

but we do not use this later. 

The following gives an algorithm to compute the weight enumerator polynomial of P 1 - 
by enumeration of P. 

Theorem 2. Let P be a subgroup o/X* . For X G X* ; let Xj G X n 6e its i-th row, namely, 
X = (xi, . . . ,x s ) T G X*. Then, the weight enumerator polynomial of P is given by 

ns s 
a=0 " ^ > XeP i=X 

where the polynomial p(h; z) (h being an integer between and n) is defined by 



p{h;z) :-- 



i + YT a =x za ( ba ~ h ) ifh = 0, 



1 + Y^Zi z a (b a - b^ 1 ) - z h b h ^ ifh>0. 

Suppose that #(P) = b m . Then, P is a strict (t,m, s)-net with 

t = m + 1 — min{a|l < a < ns, N a ^ 0}, 

where for N\ — ■ ■ ■ — N ns = we set min{a|l < a < ns, N a ^ 0} = ns + 1. 

The last condition is satisfied only if s = 1, m = n and P = X*, and then P is a 
(0, m, s)-net. 

Proof. For any function / : — > C[^], its Fourier transform / : X* — > C[z] is defined by 

K€K° n 

where K • X G T C C. Orthogonality of the characters stated in Lemma [1] implies that 
for a subgroup PcX^ 

— -r- V (K • X) = 1 if K G P x and = if Ji £ P x . 

Thus, we have the so-called Poisson summation formula 

#^yE/w= E i^E/W'-^- E /<*)■ 

" v ; XeP A'eK» " v 7 XeP ifeP x 

Now we put /(If) := z^ K ' G C[z]. Then the right most end is the weight enumerator 
polynomial of P x . On the other hand, for X = G X*, 

s I n 

fix) = ^ {k \k.x) = n Yl ^••••'^ )) n*i(£u; 

KeK» i=l V fe=(«i J ... ) /c n )6Kn J =1 



We compute the i-th component; namely, we put £j := and we prove in the following 
that the polynomial (which appeared in the above product) 



n 



g((6,..., £„);*):= E ^ U K &) 

(fti,...,K n )gK n i=i 
is equal to p(h; z), where h = fi*(xi). The coefficient of z a of q((£i, ■ ■ ■ , £ n ); z) is 

n I \ a— 1 

e n^-fe)- e «&) iie^)- 

fe=( K i,...,/t n )eK„ J'=l \kgG*-{0} / j=lKGG* 

^t(fc)=a 

This is zero if there is £j 7^ with 1 < j < a — 1, (b a — b" 1 " 1 ) if £1 = • • - = £ a = 0, and 
-fr 1 " 1 if £1 = • • • = £ a _! = 0,£ a ^ 0, because E kG g*-{o} k (0 = (E KeG * - 1 = 

6 — 1 or — 1, according to £ a = or not. Thus, <z((£i, . . • ,£ n ); depends only on the 
value . . . , £ n ). Namely, for x = (£1, . . . , £ n ) we have z) = p(fi*(x); z) where 

f l + ^ =i ^_ 6 a-l) if /» = 0, 

Pl ' Zj ' \ 1 + Eti 1 ^ a ( foa - & a ~ 1 ) - 2 h 6 fc_1 if /i > 0. 
Now for X = (xi, . . . , x s ) T G X*, 

f(x) = Hp(»*( Xi y,z) 

i=l 

and hence 

ns 1 s 

a=0 W ^ ' X=(xi,...,x s ) T eP *=1 

□ 

In the above proof, we showed: 
Lemma 4. The Fourier transform o/K n — > C[z],k i— >■ z^ fe ) zs 

= p(fi*(x), z). 

Indeed, we have a generalized version of Theorem [2J as follows. Define a function F 
from to the polynomial ring of sn variables: 

s 

F : K s n ^ C[z id \l < i < s,l < j < n], (k u ...,k s ) \-> ]\z iM . 



i=i 



A generalized weight enumerator polynomial for P 1 - is defined by 

GWpx(zy) := E 
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The Poisson summation formula then yields 

^px(^) = -1-^F(X). 



Here, by definition, 



This equals 



K=(k 1 ,...,k a ) T £K? l 



i=\ ki£K n 

Now, we define an (n + l)-variable polynomial p(/i; z , . . . , z n ) by the summation 

n 

p(h] Z ,..., Zn) := Y Z a T a(h), 
a=0 

where T a (h) are integers defined as the coefficients of z a in the polynomial p(h;z) in 
Theorem [2j Then we have the equality 

Y z m( k »x) = p(ji*(x); Zq,..., z n ). 

(The same proof as for Theorem [2] works for proving this equality. In fact, p(h; z) = 
p(/i; z°, z 1 , . . . , z n ) holds.) Thus, we proved 

Theorem 3. With the notation above we have 

1 s 

GW P x[z id ] := Y z lAkl) ■ ■ ■ z s>Kks) = ^— ^2*[[l>(tJ.*(xi);Zi t o,...,Zi in ). 

If we specialize the variables by putting Zi >a := z a in GWpx(zij), then we have The- 
orem |2j If we specialize z^ a := z a , then we have Trinker's version of the Mac Williams 
identity [27]. 

The following corollary shows that, to compute t-values of all the projections of P, 
it suffices to compute a specialization GW P ±(zi, . . . , z s ) of GW P ±(z id ) obtained by the 
substitution z id <(— z\. That is, 

GW P ,(z u ...,z s ):= Yl z^---z^ (3) 

K=(k x ,...,k s ) T 



1 S 



Corollary 1. Let P C X* . Take u C {1, . . . , s}, u ^ 0. We may consider the projection 
P u of P to the coordinates in u, which is the image of P by the projection X^ — > X^ (here 
X^ is the set of mappings from u to X ny ). Then, the weight enumerator polynomial of P u 
is obtained from GW P ± [z^j] by substituting z id <— z^ for i G u, Zi o 1 for i ^ u, and 
Zi j <— for i ^ u and 1 < j < n. Or equivalently, by substituting Zi <— z for i G u and 
Zi <— for i u in GW P x{z\, . . . , z s ) defined in (j3J). 
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Proof. We have a commutative diagram of two short exact sequences of abelian groups: 

-> P x ->■ -> P* ->■ 

t t t 

-> PJ- K£ -> P u * 0. 

The middle vertical arrow KJJ — >■ is obtained by supplementing on every i-th row for 
i (£l u. The middle and the right vertical arrows are injective. Hence, by diagram chasing, 
we have P„ = P 1 - D KJJ. By definition of GW P ± we have 

By P„ = P- 1 PlK^, the above summation is over if = (ki, . . . ,k 8 ) £ P -1 satisfying fcj = 

for every i ^ w. Compared with the definition GW P ± [z id ] := J2K={k 1> ...,k 3 ) T eP 1 - z i,K*i) " " " z s^{k 3 ), 

we notice that GW P ± [Zij] with substitutions 

f 1 for % i u, j = 0, 
hj \ for z ^ u, 1 < j < n, 

gives the GW P ±(zij). Then, the substitution Zij «— ^ for z ^ u gives the weight enumer- 
ator polynomial of P^. □ 

From the polynomial GW P ±(zi, . . . ,z s ), one can read which projection has large t- 
values. More precisely, given an s' with 1 < s' < s, we can compute the (worst) choice of 
the s' coordinates to which the projection of P have the largest exact t- value. 

For a monomial m in C[z±, . . . , z s ], we define its support C {1, 2, . . . , s} as the set of 
indices of Zi appearing in m (without counting the multiplicity). 

Corollary 2. Let GW P ±(z±, . . . ,z s ) be as above. Let Hd be the (homogeneous) degree d 
part with respect to the total degree. Suppose that Hd 0. Take any monomial rrid in Hd 
such that the cardinality Cd of its support is minimum among the monomials in Hd- For 
Hd = 0, we define Cd = s + 1 (indeed, we may define Cd as any number exceeding s). 

For a given s' with 1 < s' < s, define d' as the minimum d > 1 satisfying Cd < s'. 
Then, the projection to the support of m<j> has the exact t-value t' :— m + 1 — d'. This is 
the largest exact t-value among all the projections P u with #(it) < s', ^ u C {1, . . . , s}. 

Proof. For a given nonempty u C {l,...,s}, we consider the projection P u . Suppose 
that P^ 7^ {0}. Let ao be the minimum degree of nonconstant monomials in the weight 
enumerator polynomial of P„. Then, the exact t-value of P u is m + 1 — ao by Theorem [U 
The previous corollary implies that the degree d part of the weight enumerator poly- 
nomial of P^ is obtained from Hd by substituting Z{ z for i £ u and Z{ for % ^ u. 
Thus, ao is the minimum positive integer d such that Hd does not vanish by this substitu- 
tion. Because the coefficients of Hd are nonnegative (being a part of weight enumerator 
polynomial), the non- vanishing property is unchanged if we substitute Z{ <— z% for ieti 
and Zi ^ for % u. For a given P^, the u with the minimal cardinality satisfying 
this non-vanishing property is given as the support of m^; because if m is a monomial 
appearing in P^, then taking u as the support of m we have the non-vanishing property. 
Conversely, if we have the non-vanishing property, then there is a monomial m in Hd 
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whose support is contained in u. Consequently, for a given s' and for all u with cardi- 
nality < s', the minimum d such that Ha has the non- vanishing property with respect to 
u is given as the minimum d such that ca < s', and by choosing u to be the support of 
m d . □ 

2.2 On computing the t-value for general nets 

In Theorem [2] we showed a MacWilliams type identity for digital nets over finite abelian 
groups. We now investigate this result when P is not a subgroup of X* , but an arbitrary 
subset of X* of size b m . 

For point sets V = {x , . . . , cc^m-i} c [0, l] s it is known from Hellekalek [8] that V is 
a strict (t, m, s)-net in base b where 

t = m- max |z > : Vfc G {0, . . . , b m - 1} S with fi(k) < z : 

7— E Gwal fe (cci) = / G wal fc (x{)daj[. 
6 *=o ■ / r°-U s 



This means that strength of the net m — t equals the Walsh degree of exactness of the 
QMC rule based on the point set V . By the explanation after Definitional we may fix any 
n > m and let P C X^ be the approximation of V which is a combinatorial (t, m, s)-net. 
Then, the above formula is equivalent to 



t = m- max jz > : VJi G with /j(K) < z : 



Note that the last term is 1 for K — and for 

For a general point set P one can still compute the polynomial 

s ms 

1^ e n^^)^)=E^ 

™ ' i=(n,..,x s ) T eP «=i «=o 



where 



To see this note that 



#(p) ^ v y ' y #(p) 

7rv 7 x=(n,.,x a ) T ePi=i 7rv 7 x=(x 1 ,...,x s ) T eP 



^ / x=(x 1 ,..,x s )T eP ra» 



z^ # p 
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For digital nets we have N a = N a since the sum ^2 X€P K*X takes on only the values 
1 if K G P 1 - or otherwise. For general point sets this does not hold anymore. 

On the other hand, for general (t, m, s)-nets, N a * ^ for some a* > 1 implies that 
there is a K 6 K. s n with fi(K) = a* such that J2xgp K • X ^ 0. Thus it follows from 
(J4]) that in this case 

t > m + 1 -min{a|l < a < ns,N a ^ 0}, (5) 

where the minimum is defined as ns + 1 if Ni — ■ ■ ■ — N ns = 0. This may only occur if 
s — 1, m — n and P = X^, in which case P is a digital (0, m, l)-net. 

Thus the method from the previous section can still be used to obtain a lower bound 
on the t- value for (t, m, s)-nets in base b. 

Note that for general point sets equality in §5§ may not hold. To see this, note that 

. 6 m -l s 
1=0 i=l 

only depends on x\ through the values fi*(xi t i). Thus we can construct an example where 
the inequality is strict in the following way. Take a strict digital (t, m, s)-net in base b with 
t < m — s — 1 [U [T71 122]. Then shift all points in the elementary interval n^ =1 [l — 1) 
to the point (1 — ft -1 , . . . , 1 — leaving the remaining points unchanged. This new 
point set has a t- value t > m — s, since the interval [1 — b~ 2 , 1) x P3^ =2 [1 — ft -1 , 1) is empty. 
But the polynomial 

j b m -l s 
1=0 i=l 

stays unchanged, since the values of n*(x^i) are the same for both point sets. Thus, the 
quality parameter t' for the new point set satisfies 

t' >m — s > m — s — 1 > t — m — min{a|l < a < m, N a ^ 0}. 



2.3 The first algorithm for computing the lvalue and weight 
enumerator polynomial 

Theorem [2] yields the following algorithm for computing the exact t- value of a digital net 
(or finding a lower bound on t for general point sets). The algorithm works for digital 
nets over finite abelian groups, and consequently for finite rings. 
A straightforward computation shows the following: 

Lemma 5. For h > let p(h; z) be the polynomials given in Theorem^ For < h < n, 
we have 

and 

1 _ (h?>) m+1 
p(0; z) = (1 - z) . . + b m z rn+1 . 

1 — (bz) 



14 



Note that if one wants to compute the i-value, then it is sufficient to compute N a 



for 1 < a < m, thus we can replace p(0;z) = (1 — z) 1 j~7^ V b m z m+l by p(0; z) 

(1 ~~ z ) 1 ~i-(L) + • Then to compute the t- value, the computation of 



^ElWto);*) (6) 

" xePi=i 
can be done using the following formula 

(T^)pEnfl-w w ), 

v ^ leP i=i 

where 

/i*0i) if/i*(xi)>0, 



We have the formal expansions 
1 - z 



m + 1 if n*(xi) = 0. 



1 - bz 

and 



1 + (6 - l)z + (6 2 - &K + • • ■ + (6 m - r" 1 )^ + 



H) =e*° e n^-^-'j). 

/ a=0 C1 c>0 i=l 



ci,...,c s >0 
ciH hc s =a 



Since it is sufficient to compute ([6]) only up to degree m for computing the t- value, it 
suffices to use the polynomial 

Q m (z) := [1 + [b - l)z + (b 2 - b)z 2 + --- + {b m - b m ~ l )z m ] s (mod z m+l ) 

m s 

=$> a e n^-L^-'j)- 

a=0 ci,...,c s >0 i=l 
ciH hc s =a 

Algorithm 1. i. Given: digital net V = {xq, . . . ,Xf,™-i} C [0, l) s fewer a finite 
abelian group G, finite field or finite Frobenius ring with b elements). Let Xi = 
( x i,i, • • • j Xl,s)- 

2. Compute the coefficients N a of z a for < a < m of the polynomial 

b m -l s 

Qmiz)^ E II C 1 - Q>*f M ) (mod (7) 

1=0 i=l 

where 

\- log;, x] if x> 0, 



m + 1 if x = 0. 



5. Let the coefficient of z a of the weight enumerator polynomial (J2j) be N a . Then 
compute 

t = m — min{a|l < a < m, N a ^ 0}, 
where the minimum is defined to be m + 1 if Ni — ■ ■ ■ — N m = 0. 
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4- Return t. 

Remark 2. If one uses Algorithm^ for a general point setV = {xq, . . . , av-i} C [0, l) s ; 
then the returned value t is a lower bound on the quality parameter of the point set, i.e., 
it implies that V is not a (t — l,m, s)-net in base b. 

Remark 3. 1. Note that the proof of Theorem^ can be modified by setting 

/(*)== E f(K)(K.X), 

where one can choose £ > m — t (for instance if to is a known lower bound for the 
given net, then one can choose £ — m — to). One then obtains the polynomials 

pt(l;z) = l-z, 

Pi(2;z) = 1 + (6- l)z -bz 2 , 

Vt(t\ z) = 1 + (b - l)z + (b 2 - b)z 2 + ■■■ + (b^ 1 - b 1 - 2 ^- 1 - b l ~ l z\ 
pi(0; z) = l + (b-l)z+ (b 2 - b)z 2 + ■■■ + (b^ 1 - b £ - 2 )z e ~ l + (b e - b e ^)z e . 

and one computes the polynomial 

j b m -l s 
(=0 i=l 

where (^(x^i) = \—\og b x], and 

t = m + l- min{a| 1 < a < £ : N a ^ 0}, 

where the minimum is defined to be i + 1 if N± — ■ ■ ■ — Ng — ( since we assume 
that £ > m — t, N\ = ■ ■ ■ = Ng = can only happen if the exact t-value is m — £). 
Again, this computation can be simplified using Lemma\^ 

This way one can reduce the computational cost of calculating the weight enumerator 
polynomial modulo z m+1 . 

2. On the other hand, if one wants to compute the t-values of a point set ofb mu points 
for several values ofm£ {mi, . . . , m u } (mi < m 2 < • • • < m u ), then one can, for in- 
stance, choose £ = m u for all cases. In this case, the sum ^=o _1 Yli=iPt{f J '*{ x i,i)'i z ) 
can be reused when computing Y^i=q 1 Y\l=\ Pt{^*{ x i,i)\ z )> *.e., one only needs to 
compute Y^i=b™r _1 ili=i Vlil 1 * { x lj)i z ) an d a dd it to the previous result for the sum 

£S _1 rG=iP*(/**(3M);*)- 

We now discuss the case when computing N a for < a < D with D > m, i.e., the 
computation of the weight enumerator polynomial. Using Lemma [Q we can improve the 
order of the computational complexity of computing the weight enumerator polynomial 
as follows. Consider the special case where 

Any row of X G P \ {0} is non zero. (8) 
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Then in the computation of 

4b E iliWbM (9) 

" xeP\{o} »=i 

no p(0; z) is involved. Thus, we may compute (EJ) as follows. Let Z be a new variable, 
which will be substituted by bz later. Compute 



-1 s 



^ / 1=1 t=i 

Note that dividing a polynomial by 1 — Z is an easy task. Then 

WPpx(z) = (1 - zY^R{bz) + -Lp(0; z) s 



Now we consider the general case. In computing rii=iP(/ i *( a ';,«)! z ) m 

Z=0 i=l 

we count the number £{x{) of i with n*{xi t i) = 0. We prepare s memory for the 
polynomials. The 0th one is accumulating the sum of YLi=iP{f J, *{ x i,i) : i z ) f° r x i with 
£(xi) = 0, by the method described as above. The rth one is accumulating the sum for 
xi with £(xi) = r. The product is separated into two parts, the product of those with 
fi(xi,i) > (for which we can use the same trick based on new variable Z as above) and 
the product of those x^i with fi(xi t i) = 0. We can factor out the latter (since the number 
is r), and may sum the former terms (with fixed r). After exhausting all < / < b m , 
finally, we can substitute Z := bz for each of the s polynomials, and add them up. 

Note that the above approach for computing the t-value can also be used for the 
generalized Mac Williams. 

We now investigate the complexity of computing the t- value via the weight enumerator 
polynomial ^ YI^q 1 Yli=iP(^*( x i,i)^ z )- To do so, the coefficients of the product n!=i(l — 
Q) Z y*( x i,i)^ only need to be computed up to z m (where the coefficient is bounded by fr ms+m , 
which can be obtained by estimating WP P ±(1)) and therefore can be computed in 0(m) 
operations. Thus the t- value can be computed in 0(Ns log N) operations. 

The computational complexity for computing the minimum distance in a binary linear 
code is iVP-hard [29] in terms of the dependence on the dimension. This implies that 
computing the t-value of a digital net in base 2 is also iVP-hard. 

1 Notice that the degree of the polynomials is in practice up to, say 30 (which yields 2 30 sa 10 9 points). 
Even if one requires more detailed knowledge of the weights, the degree of the polynomials is likely below 
several hundreds. Thus polynomial multiplication algorithms by Schonhage and Strassen |23j . Cantor and 
Kaltofen [T] or the Toom-Cook algorithm are probably not beneficial for these computations, although 
asymptotically they have a better performance. 
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2.4 Using an inverse Mac Williams identity to compute the t- 
value 



Although we stated the MacWilliams identity for digital nets and the dual group, it can 
be understood as a relationship between a group and its dual group. As such, it is quite 
obvious that a MacWilliams identity as stated in Theorem [2] in the reverse direction 
is possible. However, the relationship of the weight enumerator polynomial with the 
t- value is lost. In the following we prove that for a certain choice of weight on the 
other hand, there does exist a relationship between the sum of products of polynomials 
(&sr Yln^ 1 rii=i P(^*( x n,i)', z) in Theorem [2]) and the t- value for the reverse statement of 
the MacWilliams identity. This is shown in the following theorem. 

Theorem 4. Let G be a finite abelian group consisting ofb elements. Let P C X* be a sub 
group of size b m . For X G X^ ; let Xj G X n be its i-th row, namely, X = (xi, . . . , x s ) T G X* . 
Then V is a strict digital (t,m, s)-net over G with 

t =(1 - s)(m + 1) + deg ^ - (1 + {b - l)z + • • • + {b m - b m ~ l )z m - b m z m+1 ) s 

s 

ysm-rn ^ ^(m+l)s _ -Q 

X&P i=l 

where v*{xi) = for Xi ^ and v*(Q) — m + 1. 

Proof. For k G K m let 

q(k] z) = z m+1 - u * {x) (k mx)- z m+l J2 k * x - 

m 



If k — we have 

m 

p(0; z) := g(0; z) = z m+1 ~ u ^ x) - z m+1 ^1 = 1 + ^(6 Q - b a ~ v )z a - b m z r 



xGX m xeX m a=l 

Assume now that k ^ 0. Then 



g(fc;z) = 1 + J2 Z ° Yl k»x- z m+1 k * x 



a=m+l — i/* (&') 



Since fc ^ we have V* fc • x — 0. We consider now the double sum. Let k 
. . . , « m ) and cc = (£i, . . . , The condition a = m + 1 — v*(x) implies that x 
(0, . . . , 0, € m +i-a, ■ ■ ■ , £m) with Cm+i-a ^ °- Then 



a-1 



^ ] k»X— I ^ ] K rn _|_i_ a (^ m +i__ a ) J J ^ ] K m+l-i((m+l- 

XSX m \5m + l-aGG-{0} / 1=1 Cm+l-iGG 

a=m+l— 

if /tj 7^ for an i G {m + 2 — a, ... , m}, 
6 a - fr 1 " 1 if n m = ■ ■ ■ = K m+l _ a = 0, 

—b a if K m = ■ ■ ■ = K m+ 2- a = 0, K m+ i_ a 7^ 0. 
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Thus g(/c; z) depends only on and we have 

p(ji(k);z):=q(k;z) = l+ £ (6° - & a ~V - 6 

a=l 

Thus, for < /i < m we have 

deg(p(/i; z)) = m + 1 — /i 

and the coefficient of z m+l ~ h of p(h; z) is negative in all cases. 
Let 

s 

Then, since the leading coefficients of the polynomials z) are always negative, it 

follows that 



deg(Q) = max deg 1 xpijiiki)] z) 1 = s(m + 1) — min n{K), 
which is the same as 



min u(K) = s(m + 1) — deg(Q). 

-fsrGP x -{o} 

Thus Theorem [T] implies that 

£ = m + 1 — min (J>(K) = (1 — s)(m + 1) + deg(Q). 

A"gP x -{0} 

We now find a quick way of computing the polynomial Q. We have 



xeP i=i 

Thus 



Q(*) = _(p( ; z)) s + #(P-L) ^ z( m+1 ) s [J(Z-^^ - 1). 



xeP i=i 



□ 



2.5 The second algorithm for computing the lvalue 

We now present an algorithm for computing the t-value of a digital net based on Theo- 
rem HI Again the algorithm works for digital nets over finite abelian groups, and conse- 
quently for finite rings. 
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Algorithm 2. 1. Given: digital net V = {xq, . . . , Xf,m-i} C [0, l) s (over a finite 
abelian group G, finite field or finite Frobenius ring with b elements). Let X\ = 

. . . ,Xl >s ). 

2. Compute the coefficients of z a for (s — l)(m+l) < a < s(m + l) — 1 of the polynomial 

b m -l s 

Q{z) = -(l+(6-l)z+- • .+(b m -b m - 1 )z m -b m z m+1 ) a +b am - m \\{z^ bm) -z m+l ), 

1=0 i=l 

where /i(0) = and for a positive integer u > 1 we have fi(u) = 1 + [log b u\ . 

3. Compute t = (1 - s)(m + 1) + deg(Q). 
4- Return t. 

As opposed to the first algorithm, Algorithm [2J is not extensible as explained in Re- 
mark [3J item 2). The value of m can however be adjusted in a similar manner as in 
Remark [3] for the first algorithm. If doing so one needs to use n(xijb m ), where m! is the 
new depth. 

The computationally most expensive step in the algorithm is 2., where products of 
polynomials need to be computed for < / < b m . Computing the product in general 
requires 0(s 2 m) operations. However, only the m+1 most significant coefficients need to 
be computed to obtain the t- value and these m + 1 coefficients can be computed in 0(ms) 
operations. This can, for instance, be done by computing the first m + 1 coefficients of 
the product of reciprocals 

s 

i=l 

which are the same as the m + 1 most significant coefficients of Y[i=i( zlJ, ^ Xl ' ibm ^ ~~ z m+1 ), 
and then storing them in the correct place. Thus the computational cost of the algorithm 
is O(NslogN) operations. 

In Section 14.31 we will discuss the relation between our character-theoretic dual P 1 - and 
the dual net defined by a ring-theoretic inner product introduced in [16]. It is proved 
that for a wide class of finite rings these notions coincide in a suitable sense. However, 
there are counter examples for general finite commutative rings R and P in which the 
ring-theoretic dual is strictly larger than the character-theoretic dual. On the other hand, 
it is proved that if P is a free -R-module, the formula in Theorem [T] holds if P 1 - is replaced 
by the ring-theoretic dual net. 

3 Numerical result 

As a proof of concept we computed the i-values of digital nets obtained from a Sobol' 
sequence as implemented in Matlab 2011a. The results are presented in the table below. 
In these experiments, Algorithm [2] is slightly faster than Algorithm [1] (although both 
algorithms have not been optimized in our experiments). 

The table also contains the exact t-values of the Sobol' sequence as computed in [2] for 
dimension up to 10 in the first column. We remark that using different direction numbers 
changes the t- values for particular m and s, see [2]. 
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Table 1: The exact t- values for digital nets obtained from Sobol' sequence as implemented 
in Matlab 2011a. The last row (marked by ts) contains the t- values of the Sobol' sequence 
from [2]. 
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4 Nets over finite abelian groups and over finite rings 



In this section we generalize the definition of (t, m, s)-nets and study digital nets over 
finite abelian groups and show how this theory relates to digital nets defined over finite 
rings. 

4.1 A generalization to (T,M,s)-nets 

As considered in Section 121 our framework is a point set P C X* = M s n (G), which is 
injectively mapped by ip" 1 : X* — > [0, l) s into the s-dimensional cube. 

Note that the cardinality of a subgroup P of X* is not necessarily of the form of b m , 
and the cardinality of the points in the elementary 6-adic intervals need not be of the 
form of 6*. This leads us to the following definition. 

Definition 7. Let M be a positive integer and for any integers di, . . . ,d s > 0, let Y := 
Y (di, . . . , d s ) denote the union of 

{(l,l),(l,2),...,(l,d 1 )},{(2,l),(l,2),...,(l,d 2 )},...,{( S ,l),(l,2),...,(l,4)}. 

A finite subset P C X* of cardinality M is a (combinatorial) (T, M, s)-net in base b 
if and only if for any b dl+d-2 ^ hds < M/T, the composition P C X* — > X* (Y) is uniform. 
IfT is the smallest value such that P is a (T, M, s)-net in base b, then we call P a strict 
(T, M, s)-net in base b. 

In this case, for any integer d < M/T, the number of points in each elementary b-adic 
interval of volume b~ d has the same number (depending only on d, namely M/b d ) of points 
in <^ _1 (P). 

In this notation, the original (t,m, s)-net is a (&*, b m , s)-net. On the other hand, let 
P be a (t, m, s)-net in base b for some admissible parameters t, m, s, b. Let A > 1 be an 
integer and let the multiset Pa be the point set where each point is taken with multiplicity 
A. Then P\ is a (A6*, Xb m , s)-net in base b. If P is a strict (t, m, s)-net in base b, then P\ 
is a strict (A6*, Xb m , s)-net in base b. For A not of the form b k , this provides examples of 
point sets satisfying Definition [71 

4.2 Generating basis and heterogenous case 

Let P be a subgroup of X* . In practice, we need to enumerate the points in P. A possible 
way is to find a generating set of P as an abelian group. More precisely, the structure 
theorem of a finite abelian group states that 

r 

i=l 

where qi is a prime power. Then, we precompute the r elements of P corresponding 
to (1, 0, . . . , 0), . . ., (0, 0, ... , 0, 1). By a lexicographic enumeration of the elements of 
f|[ =1 we can enumerate points in P through the above isomorphism. In practice, it 
may be desirable that for any 1 < r' < r, the image of Y[ i= i ^qt has a low discrepancy 
property. Then, we can apply QMC-integration for the first Yll=i ^ points, and if the 
result is not satisfactory, we may increase r', similarly as for (t, s)-sequences [H Chapter 4]. 
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We remark that in the results so far, we may choose different (i.e. non-isomorphic as 
abelian groups) G of the same cardinality b in the coordinate of X n = M s>n (G). Namely, 
we may choose n x s possibly different finite abelian groups Grjj. By taking its dual G*j, 
all the results so far hold. 

We may even change the cardinality of Gij, so that bij := #(Gi j) depends on 
The embedding of the z-th row YI™ =1 into the z-th coordinate in the interval [0, 1) 
is given as follows: first equating [0, 1) into b^i intervals, then each interval is equated 
into bi t 2 intervals, and so on. For a subgroup P C rTi<i< s i<j<n we can define 
P 1 - C Yli<i< s x<j<n^i,j- ^ e can define elementary (frjj)-adic intervals. For a subset 
Y C {1, 2, . . . , s} x {1,2,..., n}, we define its co-volume Vol(F) := Ylu 7 -) e y 

A point set P is then a (T, M, s)-net if the cardinality of P is M and the mapping 
from P to Ylu j) € y is uniform for any Y of the form Y(d\, . . . , d s ) whose co-volume is 
smaller than or equal to M/T. 

A version of an NRT-weight and a MacWilliams-type identity can be defined and 
proved, but we omit its explicit description. 

4.3 Ring-theoretic dual net 

The original notion of the dual net [IB] is defined in the case where the digital net is 
defined using a finite field, and the definition uses an inner product. We introduce a 
straight forward generalization to a finite ring. We use the letter R for finite rings in the 
following and consider digital nets over R (rather than G). Hence in the following we use 

X n :={x = (^,...,C n ):i t eR}. 

Definition 8. Let R be a finite ring. Then, X^ is a free left R-module of rank sn. Let 
P C X^ be a sub R-module. If P is a (t,m, s)-net in the sense of the combinatorial 
(t,m, s)-net (Definition^), then we call P a digital (t,m, s)-net over the ring R. 

This definition is slightly more general than that given in [10], which treats the case 
of a free R- module P. 

Recall that P C X* is free of rank m if there are X\ , X 2 , . . . , X m e P such that every 
element X G P is uniquely represented by an (x±, . . . , x m ) G R m as 

X = x±Xi + • • • + x m X m . 

Namely, Xi, . . . , X m are linearly independent over R and generate P as an .R-module. 

Definition 9. Let C±, . . . ,C S G M n ^ m (R) be matrices and let Xi, . . . , X m be defined by: 
the j-th row of Xi (n- dimensional) is the transpose of the i-th column ofCj. Assume that 
Xi, . . . ,X m are a free basis. Then we call the set PcXJ which is generated by the free 
basis Xi, . . . , X m , the digital net generated by the matrices C±, . . . , C s . 

The above P is the same point set as defined in [10] in Section 1. This follows by just 
comparing the columns of Cj with the rows of X m . The assumption that #(-P) = b m is 
equivalent to that Xi, . . . , X m give a basis of P. 

When R is a finite field, the notion of the dual net is introduced in [TH]. A straight 
forward generalization to a finite ring is: 
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Definition 10. Let P C X* be a left R-module and let 1,1' 6 P with X = (xi t j) and 
X' = (x'^). We define 

p v : = {x' g X s n \J2 x 'ij x ij = f° r an v x eP}a x*. 

hi 

The expression of the exact t- value of P in terms of the minimum NRT-weight of P v 
is given in [TO] in the case where R is a finite field. 

To deduce their result from our character theoretic dual P 1 - G X^, we need one 
definition on finite commutative rings, cf. |31j . 

Definition 11. A finite commutative ring R has a generating character if there exists a 
character c : R — > T as an additive group such that the pairing obtained by the composition 
with the multiplication fi 

R x R A- R A- T 

is a perfect pairing, namely, 

9 C : R^ R*,x ^ c(x ■ (-)) 

is an isomorphism of abelian groups. We call such a ring R a finite commutative ring 
with generating character. 

In [31] it was shown that finite commutative rings with a generating character are 
equivalent to finite commutative Frobenius rings, which in turn are equivalent to finite 
commutative Gorenstein rings. 

It is easy to see that the class of such rings is closed under taking cartesian products 
and subrings (because of the finiteness: a subring S C R has S — > S* as a restriction of 
R R*, which is an injection and hence bijective). It includes Z&, since we may then 
take c : Z5 — > T as an injective group homomorphism. It includes finite fields, since a 
finite field R is a vector space over a prime field Z p , we may take any nontrivial Z p -linear 
mapping R — > Z p and an injective group homomorphism Z p — > T, and then c is obtained 
as the composition c : R — > Z p — > T. Thus, it is a large class of finite commutative 
associative rings, but it does not include all finite rings, see Remark |4] below. 

Lemma 6. Suppose that R is a finite commutative associative ring with generating char- 
acter. Then the isomorphism 

of abelian groups induced by 9 C componentwise, induces an isomorphism 

P v -> P x . 

Proof. For X' e X s n , 6 C (X') = {0 c {x' id ■(-))) G K*. The condition X G P y is equivalent 
to 

x'i jXij = for any X G P. 

hi 

The condition X G P 1 - is, through the identification R — > R* above, 

c C^ x 'i,j x i,j) = 1 f° r an y x e p. 

hi 
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Thus, P v C P is automatic. For the converse inclusion, we need the fact that P is an 
P-module. The set 



I hi 



is then an P-module (by commutativity of R). The condition X' e P -1 is equivalent to 
c(S') = 1. This implies that any element s E S has the property 6 c (rs) = 1 for all r. By 
the assumption of the character, this implies s = 0, hence S = 0, and X' £ P v . □ 

Corollary 3. Let R be a finite commutative ring with generating character. The result 
in ITby based on the ring theoretic dual net P v follows from our character theoretic result 
Theorem [3 

Remark 4. Let R be the finite commutative associative ring 1* p [x,y]/(x 2 ,xy,y 2 ). Then 
the Itp-linear span P :=< x >C R is an R-module of cardinality p and satisfies P v =< 
x,y >, whose cardinality is p 2 . On the other hand, P ± has the same cardinality p as P. 
Thus, R does not have a generating character and the cardinality of P v is strictly larger 
than the cardinality of P -1 . 

We see that if R is finite commutative associative ring with generating character, then 
the two notions of duality coincide, and the minimum NRT-weight of P v gives the strict 
t- value. 

Suppose that R is a finite commutative associative ring, not assuming that it has a 
generating character. Still, if P is a free P-module, the same result for the minimum 
NRT-weight holds, by the following proposition. 

Proposition 1. Under the above assumption on R, let P C X* fe a free R-module. Then, 
for any Y , the composition P C X* — > ~K s n (Y) is surjective if and only if P v fl X* (Y) = 0, 
where X s n (Y) := {X = (xij)\x id = for all £ Y} C X s n . 

Proof. For an P-module V, let us write P := Hom^fV, P). It is obvious that if / : R n — > 
R m is surjective, then / : R m — > R n is injective. The converse, / injective does not 
necessarily imply / surjective for a general ring, but it is true for a finite ring P (see the 
next lemma). 

This makes it possible to replace P* with P, then P 1 - with P v , in the proof of Lemma[3l 
if P is a free P-module. □ 

Corollary 4. Suppose that P is a free R-module. Then, the statement of Theorem^ is 
true if P 1 - is replaced with P v . 

The following result is a deeper result than it looks like. 

Lemma 7. Let R be a commutative associative ring of Krull dimension zero ( this con- 
dition is satisfied if R is a finite commutative associative ring). Then, if a morphism 
f : R n — > R m of the R-module is injective, then f is surjective. 

Proof. For such a ring, it is proved in [131 Proposition 2.6] that for any mxn matrix A of 
coefficients in P, if the multiplication of A from left to R n is injective, then one can extend 
A to an invertible matrix by adding m — n columns. In particular, the multiplication of A 
from the right to a horizontal vector R n is surjective. This means that the multiplication 
of the transpose T A from the left is surjective, which is the representation matrix of /, 
hence / is surjective. □ 
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